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Abstract 

Wc present a distance measurement for the semiregular variable S Crateris (S Crt) based on its annual 
parallax. With the unique dual beam system of the VLBI Exploration for Radio Astrometry (VERA) 
telescopes, we measured the absolute proper motion of a water maser spot associated with S Crt, referred 
to the quasar J1147— 0724 located at an angular separation of 1.23°. In observations spanning nearly 
two years, we have detected the maser spot at the LSR velocity of 34.7kms _1 , for which we measured 
the annual parallax of 2.33±0.13mas corresponding to a distance of 430^ 25 jpc. This measurement has an 
accuracy one order of magnitude better than the parallax measurements of HIPPARCOS. The angular 
distribution and three-dimensional velocity field of maser spots indicate a bipolar outflow with the flow 
axis along northeast-southwest direction. Using the distance and photospheric temperature, we estimate 
the stellar radius of S Crt and compare it with those of Mira variables. 
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1. Introduction 

Very Long Baseline Intcrfcromctry (VLBI) is a pow- 
erful technique for obtaining positions of celestial ob- 
jects with milliarcsecond (mas) level accuracy. The VLBI 
Exploration of Radio Astrometry (VERA) telescopes are 
a Japanese VLBI array dedicated to phase referencing 
VLBI (Kobayashi et al. 2003). VERA consists of four 
20 m diameter antennas at Mizusawa, Ogasawara, Iriki 
and Ishigaki-jima (see figure 1 of Petrov et al. 2007). To 



overcome phase fluctuations and the limited integration 
time in conventional fast-switching VLBI, VERA has a 
dual beam system that allows simultaneous observations 
of the target and reference sources separated by 0.3 to 
2.2 degrees. This advanced capability of VERA can be 
used to measure annual parallaxes and proper motions of 
masers with a 10 micro arcsecond (/zas) level accuracy. 

An annual parallax gives a simple geometrical measure 
of the distance, free from the complex assumptions re- 
quired in other distance estimators. Studies of Asymptotic 
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Giant Branch (AGB) stars by Whitelock et al. (2008), 
Kurayama et al. (2005), and Vlemmings et al. (2003) are 
based on the annual parallax derived from VLBI obser- 
vations, and demonstrated the the ability of VLBI as- 
trometry. More recently, successful measurements of dis- 
tances and proper motions with VERA have been re- 
ported (Hirota et al. 2007; Honma et al. 2007; Imai et al. 
2007; Hirota et al. 2008; Sato et al. 2007). In the present 
paper, we report the observations of S Crt, concentrating 
on the distance measurement. S Crt (IRAS11501— 0719, 
IRC- 10259, AFGL4830S) is an AGB star, with a pulsa- 
tion period of 155 days (Benson et al. 1990). S Crt has a 
variability type of SRb in the General Catalog of Variable 
Stars (Antipin et al. 2005). The apparent magnitude at K 
band is 0.786±0.314mag (Cutri et al. 2003). The optical 
light curve in V band is presented in figure 1(a) exhibiting 
~0.8mag variation. 

The distance to S Crt has been measured in various 
ways. Bowers & Johnston (1994) and Patel et al. (1992) 
determined the distance to be 420 pc and 285 pc, respec- 
tively, using the period-luminosity (PL) relation of Mira 
variables. Perryman et al. (1997) used HIPPARCOS data 
to measure a parallax of 2.04±1.31 mas, corresponding 
to a distance of 490 pc, but with an large uncertainty. 
Recently, the new HIPPARCOS catalog (van Leeuwen 
2007) gives the parallax of 1.27±0.92mas, corresponding 
to a distance of 787 pc. VERA offers the potential to de- 
termine the parallax with a better accuracy so that the 
inconsistencies in these distance measurements can be re- 
solved. 

The observations and data reduction are described in 
section 2. In section 3, we present absolute positions and 
internal motions of maser spots. The annual parallax of 
S Crt is presented and converted to the absolute distance 
in this section. Finally, in section 4, we examine the prop- 
erties of S Crt in the context of the new distance measure- 
ment. 

2. Observations and data reduction 

2.1. Observations 

We have conducted a series of monthly VLBI obser- 
vations from October 2005 to May 2007 with VERA. 
The duration of each observation was typically 8.5 hours, 
yielding a net integration time of 5 to 6 hours. The 
typical synthesized beam size (FWHM) was 1.5 mas x 
0.8 mas with a position angle (P.A.) of —30°. To ob- 
tain the positions of maser spots in S Crt, the quasar 
J1147— 0724 was simultaneously observed as the position 
reference. The J2000 a priori coordinates of two sources 
are (a, 8) = (11 h 52 m 45.0981s, -07° 35' 48.072") for 
S Crt, and (11 h 47m 51.554035s, -07° 24' 41.14109") 
for J1147-0724. J1147-0724 is classified as a "candi- 
date" source in the catalog of the International Celestial 
Reference Frame (ICRF), and its position errors are 
270 /^as and 290 /xas in R.A. and Dec, respectively (Ma 
et al. 1998). The separation and P.A. of these sources are 
1.23° and -80°. 

From a total of 17 observations, we used 12 observations 



in the present study. The observation in 2005 November 
was not used because the Ishigaki-jima antenna did not 
participate. Bad weather conditions in four observations 
held between August and October in 2006 resulted in 
much higher system noise temperatures (>500K) at some 
of the stations. The shorter integration times and un- 
solved phase fluctuations in these four observations de- 
graded the image qualities of S Crt in both self-calibration 
and phase reference analyses. 

The observation status is summarized in table 1, which 
is organized as follows: Column (1) — The ID number 
of epoch. The numbers with an asterisk "*" indicates 
that the observation was used in the estimation of the 
annual parallax. Column (2) — Date of the observation. 
Column (3) — The year and day of the year (DOY) . In the 
observations in 2007/013, 2007/096, and 2007/130, the 
phase referenced images are scattered into several com- 
ponents, while the self-calibrated images resulted in a 
single bright component, clearly indicating that the at- 
mospheric conditions on these days prohibited succesful 
22 GHz phase referencing. We did not use these three ob- 
servations for the parallax estimation. In parallel with the 
VLBI program, we have monitored S Crt with single-dish 
observations at Iriki station since 2003 September with 
the typical interval of one month (see also Shintani et al. 
2008). 

A data recording rate of 1024 Mbps was adopted with 
the VERA DIR2000 recording system, which yields the 
total receivable bandwidth of 256 MHz with 2-bit digiti- 
zation. The 256 MHz bandwidth data of left-hand circular 
polarization were divided into 16 IF channels of 16 MHz 
band width, one of which was used to receive the maser 
emission and the others were used to receive the contin- 
uum emission from Jl 147— 0724. Cross-correlation was 
carried out with the Mitaka FX correlator (Shibata et 
al. 1998) at the National Astronomical Observatory of 
Japan (NAOJ). In most observations, the IF channel as- 
signed to the water maser was divided into 512 spectral 
channels, yielding a frequency spacing of 31.25 kHz cor- 
responding to a velocity resolution of 0.42kms~ 1 . The 
observations in 2006/069 and 2006/129 provide two times 
higher frequency resolution than other observations, yield- 
ing a velocity resolution of 0.21 km s -1 . For the data of 
J1147— 0724, the IF channels were divided into 64 spectral 
channels in all observations. 

2.2. Calibration and Imaging 

We used the Astronomical Imaging Package Software 
(AIPS) developed in the National Radio Astronomical 
Observatory in the data reduction. Amplitude calibra- 
tion was achieved using the system noise temperatures 
and gains logged during the observations at each sta- 
tion. In the fringe search process of the reference source 
J1147— 0724, we used the task fring with a typical inte- 
gration time of 2 minutes, with solutions of fringe phases, 
group delays, and delay rates obtained every 30 seconds. 
Using the task TACOP, these solutions were transfered 
to the data of S Crt in order to calibrate the visibility 
data. Phase and amplitude solutions obtained from self- 
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Fig. 1. Light curves of S Crt. (a) The optical light curve at V band provided by the ASAS (Pojmanski 1997). (b) The water mascr 
light curve at 22 GHz obtained with the single-dish monitoring program with VERA (Shintani et al. 2008). Peak flux densities 
of major velocity components are presented with different symbols. The symbols of "+", "O", "x", and "□" indicate the 

peak fluxes of the individual components at the LSR velocities of 34.7, 36.9, 38.5, 39.5, and 40.0kms~ , respectively. There were 
no available data before the year of 2003.8. The hatched area indicates the period of the present VLBI monitoring observations. 



calibration of the Jl 147— 0724 were also transferred to the 
S Crt data. 

Since the delay-tracking models used to estimate a pri- 
ori delays in the Mitaka FX correlator were not accurate 
enough for astrometry with VERA, we have applied bet- 
ter estimates calculated with the CALC3/MSOLV soft- 
ware package (Jikc ct al. 2005). Manabe et al. (1991) 
gave a brief description of this package in the report to 
the International Earth Rotation Service (IERS). CALC3 
has slight differences in the physical models than those 
adopted in CALC which was developed by the research 
group at NASA Goddard Space Flight Center (GSFC). 
Comparing the two a priori values from Mitaka FX cor- 
relator and the CALC3/MSOLV package, we obtained the 
difference between the two estimates, then applied them 
to the visibility data using the tasks TBIN and CLCAL in 
AIPS. We note that the zenith excess path lengths due 
to the wet atmosphere measured by the global position- 
ing system (GPS) at each station are considered in the 
CALC3/MSOLV package. 

We estimated the residual atmospheric zenith delay off- 
set using the method described in Honma et al. (2007). 
This correction was applied to one of the four stations of 
VERA where the coherence of the phase-referenced image 
was maximized. Typically, the offsets were found within 
±3 cm. 



The instrumental delay caused by the difference be- 
tween two signal paths was estimated using an artificial 
noise source (Kawaguchi ct al. 2000). These delays were 
then loaded into AIPS and applied in the same manner 
as the delay tracking model correction. For each observa- 
tion, we fitted a two dimensional Gaussian model to the 
brightness distribution to find the position of the maser 
spot. Then we used these to derive the annual parallax 
and linear proper motions (/J-x, My)- 

In the single-beam VLBI imaging of S Crt, group delays 
solved with 3C279 were applied, then fringe-fitting was 
done using the brightest emission showing the LSR ve- 
locity of 34.7 km s -1 . Thus relative positions of all mascr 
spots in self-calibrated images were determined with re- 
spect to this reference spot. For each velocity channel, po- 
sition of the spot is defined as the brightness peak of the 
image. To estimate the relative motion of each spot, a lin- 
ear least-squares analysis was applied to the spot seen at 
the same velocity channel during at least two continuous 
observations. The signal to noise ratio of 10 was adopted 
as the detection criterion in the self-calibrated images. 
The rms noise of the self-calibrated images was typically 
90mJybeam _1 . On the other hand, the rms noise of the 
phase-referenced images was typically 700 mJy beam -1 , 
about one order of magnitude larger than that of self- 
calibrated images. Due to limited phase coherence the 
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Table 1. Observations 
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9110^ Drtnhpr 1 Q 
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O* 
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ZUUu / OO t 


o * 
O 


ZUUU O cLIlUdjI y o 


9nnfi /nn^ 

zuuu / uuo 


A* 
-t 




ZUUU / yi^A 


i]* 


2006 March 10 


2006/069 


6* 


2006 May 9 


2006/129 


7* 


2006 November 13 


2006/317 


8* 


2006 December 11 


2006/345 


9 


2007 January 13 


2007/013 


10* 


2007 February 21 


2007/052 


11 


2007 April 06 


2007/096 


12 


2007 May 10 


2007/130 



phase referenced images have limited signal-to-noise (30 
vs 307 in the self-calibrated images). 

3. Results 

3.1. Annual parallax and distance 

In figure 2. cross-power spectra of S Crt from 
2005/019 (upper) and 2006/069 (lower) on the Mizusawa- 
Iriki baseline are presented. The time variation of the flux 
densities with the total power spectrum are presented in 
figure 1(b) with an indication of each LSR velocity. The 
34.7kms _1 spot underwent a radio flare (e.g. Shimoikura 
ct al. 2005). This flare started in February 2006 and 
reached a maximum of 371Jy in March 31 (2006/090), 
then it decreased to 60 Jy in July. During the flare, this 
spot did not show structural change. The cross-power 
spectrum at the same time showed significantly weaker 
(~30%) flux density than that of the total power spec- 
trum, indicating that some of the flared emission is re- 
solved on the VERA baselines. 

In figure 3, we present the positions of maser spot at the 
LSR velocity of 34.7kms _1 relative to the phase track- 
ing center. Throughout the present VLBI observations, 
this maser spot was bright enough to be detected on all 
baselines in the phase referencing analyses. The proper 
motion was clearly modulated by a parallax. Based on a 
least-squares fitting analysis, the parallax was determined 
to be 2.33±0.13mas which corresponds to a distance of 
430I23 P c - Here, we adopted position errors of each mea- 
surement that are obtained as the root sum squares of 
three error factors, and the details are given in section 4.1. 

In the estimation of the parallax, we adopt a very 
small number of assumptions: the maser spot is mov- 
ing on a linear trajectory with respect to the star, i.e., 
there is no acceleration, and the reference source is fixed 
on the sky, i.e., no motions due to core shift or jet fea- 
tures. From the fitting results, the linear proper motions 
of the reference spot {nx, My) were obtained to be {^x, 
Hy) = (-1-56 ± 0.22masyr-\ -5.16 ± 0.22masyr" 1 ). 
This motion is the combination of the proper motion 
of S Crt system and internal motion of the maser spot 
in the system. Considering the offset of the refer- 
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Fig. 2. Cross-power spectra of the S Crt water mascrs 
observed with the Mizusawa-Iriki baseline on 2005/292 
(top) and 2006/069 (bottom). The stellar velocity of 
37.85 km s -1 is indicated with the vertical dashed line 
in the spectra. The blue-shifted component with re- 
spect to the stellar velocity has been brighter than 
the red-shifted one in the majority of our observations. 

ence spot from the phase tracking center (Aa,A(S) = 
(-8.49 mas, -23.81 mas), the J2000.0 absolute coordi- 
nates of this spot in 2005/292 were obtained to be (a, <5) 
= (lib. 52m 44.96969s, -07° 35' 48.0958"). This is the 
position as referenced to the position of J1147— 0742 and 
in relative offset to the original phase tracking centre. The 
uncertainty of this position is estimated as ~400 /xas based 
on the errors in our phase referencing analysis (detailed 
in section 4.1) and ICRF position of J1147-0724. 

The reference source J1147— 0724 exhibited an unre- 
solved structure. The correlated amplitude as a function 
of (u,v) distance is flat and the upper limit of the source 
size is 0.8 mas (FWHM), which is the minor axis of the 
synthesized beam in the present observation. In addition, 
we confirmed that the images of J1147— 0724 showed no 
distinctive change during the observations. 

3.2. Maser distribution and internal motions 

In figure 4, we present the angular distribution and 
three-dimensional velocity field of maser spots in S Crt 
covering a 60 mas x 60 mas region. The maser spot 
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Fig. 3. Absolute positions of the reference maser spot in 
S Crt. Filled circles indicate present results with error 
bars and crosses indicate predicted positions. Solid curve 
indicates the combined motion of the best-fit results of 
parallax and proper motion. The axes indicate the po- 
sition offsets with respect to the phase tracking center. 

with the LSR velocity of 34.7 km s -1 is placed at the 
map origin. At the distance of 430 pc, lmas corresponds 
to 0.43 AU, and 1 masyr corresponds to a velocity of 
2.04km s -1 . The color index in figure 4 shows the LSR ve- 
locity range from 34.0 to 41.0 km s . The blue- and red- 
shifted components are separated into the northeast and 
southwest parts of the area. The relative motion of each 
spot (y a , v y ) with respect to the reference spot are used 
to determine the average motion (v x , v y ) and, hence we 
obtained (v x , v y ) = (—1.605 masyr -1 , — 0.252masyr" 1 ). 
Then, we subtracted (v x , v v ) from (v x , v y ) to obtain the 
internal motions (V x , V y ), that are presented with arrows 
in figure 4. We successfully detected the internal motions 
of 26 maser spots. The typical transverse speed was ob- 
tained to be 2.72 masyr -1 , corresponding to 5. 56 km s -1 , 
by averaging the internal motions of all 26 spots. The pa- 
rameters of the motions are presented in table 2 in the 
increasing order of LSR velocity. 

For the reference spot, we subtracted the internal mo- 
tion from the proper motion (fix, Hy) and, thus, the 
proper motion of S Crt system was estimated to be 
(-3.17 ± 0.22masyr -1 , -5.41 ± 0.22masyr -1 ). In the 
new HIPPARCOS catalog (van Leeuwen 2007), the ab- 
solute proper motion of S Crt was (— 3.37± l.OOmasyr - , 
—4.67 ± 0.75 masyr -1 ), representing good consistency 
with our result within the errors. 

4. Discussion 

4-.1. Errors in the positions 

Since it is difficult to identify contributions to positional 
errors in individual VLBI observations (Kurayama et al. 
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2005; Hachisuka et al. 2006; Honma et al. 2007; Sato et 
al. 2007; Hirota et al. 2007; Imai et al. 2007), we applied 
errors in each observation by evaluating following three 
error sources: (1) a zenith atmospheric delay offset at each 
station, (2) station position errors, and (3) image quality. 
The best estimates of the zenith atmospheric delay offsets 
fell within the range of ±3 cm. Errors in the a priori 
station positions can result in additional phase residuals. 
Therefore, we assumed the offset of 3 cm even after the 
atmospheric offset calibration was done. 

For a given separation of zenith angles between two 
sources AZ, the difference of signal path lengths Al be- 
tween the two directions caused by the zenith atmospheric 
delay residual lo can be estimated as follows, 

Al = l AZ^-(^—), (1) 
aZ \cosz / 

where Z is the mean zenith angle of the observed sources. 
We adopted 1.23° as the approximation of AZ. During 
our observations, the elevation angles higher than 30° were 
dominant. We assumed Z = 50° as the zenith angle for 
estimating errors of the worst case, and obtained Al = 
0.12 cm. Using a root sum squares of major and minor 
axes of the synthesized beam (9b = 1.7 mas), we estimated 
a positional error of oz = 156 /ias from the following ex- 
pression, 

A/ 

o-z = e b x- , (2) 

*H 2 

where Xh 2 o is the wavelength corresponding to the rest 
frequency of water masers. Since the position accuracy of 
antennas are determined to be ~3 mm based on geodetic 
observations at S and X bands, we estimated the posi- 
tional error (<jd) attributed to the baseline errors to be 
5 /ias using the following expression, 

3 TTWTb 

a D = SinO s a x t , (3) 

where Osa = 1-23° is the separation angle of the source 
pair. In addition, we estimated the measurement errors 
(a i =8b/ SNR) which depends on the SNRs of the phase- 
referenced images. Finally, we obtained the position errors 
of each phase-referenced image by adding up quadratically 
the error factors, o~z,o~d, and 07. The error bars in figure 3 
represent the errors in each observation, and the averaged 
value of the errors is 167 //as. 

4-2. Maser Morphology 

Figure 4 shows the internal motions of the masers pots 
around S Crt found in the present observations. Since the 
maser spots show a bimodal distribution about radial ve- 
locities and positions, we extract the kinematic essentials. 
The analytic tools based on the diagonalization of the 
velocity variance-covariance matrix (VVCM) (Blocmhof 
2000) and position variance-covariance matrix are used. 

The elements of the VVCM matrix are presented as 

1 N 

a* = ^7— r £(Vk " Vi)(Yk,i - V k ), (4) 
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Table 2. Parameters of the detected maser spots 



ID 




X 


Y 


S 


v x 






a Vy 




(kms" 1 ) 


(mas) 


(mas) 


(Jy beam -1 ) 


(mas yr _1 ) 




(masyr - 1 ) 




1 


34.10 


-1.25 


0.34 


2.3 










2 


34.53 


-1.71 


-2.50 


2.7 


-0.01 


0.01 


4.90 


2.58 


3 


34.55 


0.76 


-5.39 


4.4 


4.15 


0.55 


-1.19 


0.06 


4 


34.68 


0.00 


0.00 


19.5 


1.61 




0.25 




5 


34.68 


-16.97 


-0.34 


4.6 


-0.55 


0.02 


0.69 


0.09 


G 


34.68 


-2.85 


-0.91 


8.1 










7 


34.72 


0.73 


-1.42 


1.8 


3.50 


0.29 


-1.09 


0.45 


8 


34.72 


1.43 


-2.00 


1.9 


2.10 


0.54 


0.88 


0.14 


9 


34.74 


1.64 


-10.71 


5.8 










10 


34.94 


-1.72 


-3.04 


1.6 


0.63 


0.50 


3.56 


0.37 


11 


34.95 


1.45 


-1.86 


35.4 


2.08 


0.89 


0.58 


0.06 


12 


34.97 


2.51 


-7.51 


2.3 










13 


34.97 


1.11 


-27.05 


2.5 










14 


35.10 


-0.03 


0.01 


8.4 


1.85 


0.74 


0.10 


0.04 


15 


35.10 


-2.89 


-0.92 


4.9 










16 


35.16 


3.52 


-2.50 


3.5 










17 


35.16 


2.13 


-5.20 


3.2 










18 


35.18 


2.38 


-7.37 


2.1 










19 


35.18 


1.14 


-27.01 


2.1 










20 


35.52 


1.30 


-2.21 


1.5 


2.05 


0.90 


0.68 


0.36 


21 


35.52 


-17.32 


-0.74 


1.3 


-0.66 


0.06 


0.80 


0.11 


22 


35.54 


-5.04 


-1.04 


0.9 










23 


36.45 


3.17 


-10.90 


1.3 


3.39 


0.52 


-0.98 


0.49 


24 


37.89 


-14.77 


-22.00 


0.8 


-4.37 





-8.49 





25 


38.32 


-27.32 


-26.17 


1.3 










2G 


38.73 


-27.06 


-28.97 


1.0 


-2.83 


0.15 


0.67 


0.45 


27 


38.73 


-19.16 


-21.19 


0.9 


-4.70 





2.46 





28 


38.74 


-29.53 


-23.72 


0.9 










29 


38.74 


-11.30 


-27.07 


1.0 


-0.55 


0.02 


0.69 


0.09 


30 


38.91 


-4.01 


-26.46 


1.1 


-1.60 


0.43 


-0.65 


0.29 


31 


39.15 


-3.67 


-25.96 


0.8 


-2.69 


0.30 


-0.30 


0.18 


32 


39.15 


-19.40 


-21.04 


2.2 


-2.67 


0.37 


0.95 


0.44 


33 


39.57 


-19.87 


-20.95 


3.4 


-3.26 





0.14 





34 


39.57 


-7.18 


-22.73 


2.2 










35 


39.57 


-19.02 


-18.52 


2.0 










3G 


39.59 


6.15 


-1.49 


1.8 










37 


39.61 


-4.73 


-28.23 


0.7 


1.49 





0.23 





38 


39.74 


-5.14 


-23.02 


1.0 










39 


40.00 


-18.91 


-18.58 


1.2 










40 


40.01 


-6.22 


-20.02 


4.7 










41 


40.03 


-4.81 


-28.14 


1.0 


1.59 





0.18 





42 


40.03 


-2.29 


-26.22 


2.0 










43 


40.16 


-9.75 


-30.09 


1.6 


0.68 


0.38 


-0.72 


0.06 


44 


40.42 


-20.47 


-20.21 


1.1 


-1.49 


0.23 


-3.01 


0.51 


45 


40.58 


-5.73 


-26.12 


2.6 


0.54 


0.05 


-0.26 


0.17 


4G 


41.00 


-5.03 


-31.63 


1.1 










47 


41.00 


-5.70 


-26.20 


5.4 


1.32 


1.27 


-0.81 


0.42 



Column (1) — Component ID. Column (2) — LSR velocity in kms -1 . Column (3) — offset positions in R.A. 
relative to the original phase center. Column (4) — offset positions in Dec. relative to the original phase 
center. Column (5) — brightness of the spot at the first detection in Jybeam -1 . Column (6) — best fit linear 
motion in R.A. in masyr -1 . Column (7) — standard errors of the motions in R.A. In case the number of 
successful detections were restricted to be only twice, we adopted the standard errors ax and cry for the 
spot as zero. Column (8) — best fit linear motion in Dec in masyr -1 . Column (9) — standard errors of the 
motions in Dec. 
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Fig. 4. Angular distribution and internal motion vectors of maser spots in S Crt. The spot color indicates the LSR ve- 
locity (see color index at the bottom). The stellar velocity of S Crt is 37.85kms — 1 . The arrow at the bottom-left cor- 
ner indicates a proper motion of lmasyr -1 , corresponding to 2.04kms _1 at the distance of 430pc. The shaded-circle indi- 
cates the size of the stellar photosphere. The dashed line represents a circular fit to the maser distribution, and has a ra- 
dius of 14.9 mas (see section 4.2 for details). The crosses indicate the stellar positions estimated from a different technique. 
The bold arrows Vi and V2 are the eigenvectors for the largest eigenvalue indicating the major axes of maser distribution. 



where the diagonal elements are the velocity dispersion. 
Here j and k denote three orthogonal space axes (e.g., 
R.A., Dec, and radial coordinate z), i denotes the ith 
maser spot in a collection totalling N(=26), and the bar 
indicates averaging over maser spots. 

The vector Vi (figure 4), which is a two-dimensional 
projections of the eigenvector for the largest eigenvalue, 
is the major principal axis of the VVCM. This axis of 
maximum internal velocity dispersion is the outflow axis. 
The axis of Vi lies at a P. A. of ~ 255°, and also make 
an inclination angle ~ 43° with the plane of the sky, with 
the southwest lobe directed into the page (away from the 
observer). The VVCM was initially evaluated in a coordi- 
nate system with axes R.A., Dec, and radial coordinate 
z increasing toward the source, and then it was diagonal- 
ized. 

We also considered the two-dimensional variance- 
covariance matrix to quantify the angular distribution of 
the maser spots. Here, the number of maser spots is 26, 
same as the case of VVCM. In the diagonalization of this 
matrix, one eigenvalue which is larger than the other by a 
factor of 3.1 was obtained. The corresponding eigenvector 



V2 presented in figure 4 indicate the spatial elongation of 
the maser distribution. The position angle of this axis is 
found to be ~ 208°. Although maser velocity and position 
axes show a discrepancy of ~ 47°, the directions of these 
model-independent axes indicate the bipolar outflow with 
the major axis along northeast-southwest. 

This dynamical property is not revealed by previous 
observation (Bowers & Johnston 1994). 

We estimate the stellar position from the present result 
by using the kinematic information. Since the motion of 
each maser spot are powered by the forming star, its posi- 
tion can be estimated as an origin of these motions. Here, 
we take one of the simplest kinematic models, which as- 
sumes a linear motion with a constant velocities from a 
single origin. This technique is essentially same as that 
used in Imai et al. (2000) and Honma ct al. (2007). Based 
on this modeling, we obtained the origin of the maser 
motions which most effectively explain the observed posi- 
tions and velocities. The stellar position from this kine- 
matic analysis, the kinematic center was obtained to be 
(X, Y) = (-7.9il.8mas, -14.9il3.3mas), which is indi- 
cated with the cross with error bars in figure 4. We used 
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the first order moment with respect to the determined 
stellar position in order to find a typical distance of maser 
spots. The obtained distance was 14.2±3.5mas, and pre- 
sented with dashed line in figure 4. At the distance of 
430 pc, the corresponding linear length is 6.11 AU. Since 
we estimated the inclination angle of the outflow axis, this 
projected length can be converted to the linear radius of 
8.96 AU. 

Nineteen years ago, the angular distribution of the 
masers was firstly observed with the NRAO Very 
Large Array in its A configuration (Bowers & Johnston 
1994). The masers were uniformly distributed over a 
40 mas x 40 mas region with the LSR velocities ranging 
evenly from 33.4 to 44.2kms~ 1 (figure 13 of Bowers & 
Johnston 1994). In their study, the blue-(red-)shifted 
spots were found in the northeast (southwest) of the maser 
region. In addition, masers with intermediate velocities 
were found in the central region. Although the extent of 
the distribution is consistent with our result, we did not 
detect the masers in central region. In figure 4, the maser 
spots are mainly found at the peripheral region of the dis- 
tribution. During the period of 2006.3 - 2007.0, we could 
not find the emission with intermediate velocities (see fig- 
ure 2(b)), and it can explain the the absence of interme- 
diate velocity maser spots . On the other hand, at the 
beginning and end of our VLBI observations, the masers 
with intermediate velocities arc found in total power spec- 
trum, and they distributed in the peripheral region. As 
Bowers & Johnston (1994) did not precisely mention the 
flux densities and sizes of each spot, it is difficult to spec- 
ify the reason for our non-detection of the masers in cen- 
tral region. However, we can speculate as to the rea- 
son; (1) Our image sensitivities were not high enough to 
detect them. The rms noise 38mJybeam _1 of Bowers 
& Johnston (1994) was lower than that of the present 
VLBI observations, 90mJybeam _1 . (2) The masers in 
central region have faded before the present VLBI obser- 
vations. (3) The masers were extended and diffuse, there- 
fore they were resolved out with the synthesized beam 
of VERA. For example, we find the effect of resolve-out 
(e.g. <30 % in 2006 March) from the comparison between 
cross-power spectrum (figure 2) and total-power spectrum 
(figure 1(b)). 

4-3. Stellar diameter and temperature 

We extract the stellar properties of S Crt using the dis- 
tance. Aringcr ct al. (1999) determined the temperature 
of photosphere (T BB ) in S Crt to be 3097±100K by fitting 
two blackbodies to its infrared spectrum. Using this tem- 
perature, we estimated an acceptable temperature range 
of 2600 - 3500 K. Here we assume that the light varia- 
tion of 0.8 mag (figure 1(a)) is same in the infrared and 
caused only by the temperature variation. Aringer et al. 
(1999) reported that Mira variables show T BB = 2418 - 
2902 K and SRVs shows T BB = 3005 - 3500 K. The esti- 
mated temperature range of S Crt is consistent with their 
finding. 

It is worthwhile estimating the stellar radius of S Crt, 
since there is no direct measurement. Using the apparent 



magnitudes in infrared (J, H, and K), T BB (3097K) and 
the distance (430 pc), we estimated the radius of S Crt 
photosphere (i?*) to be 1.81±0.14 x 10 13 cm (260±20i? s ). 
This is presented with a shadcd-circlc in figure 4, and the 
center of circle is placed on the kinematic center. Using 
this radius, we estimated an acceptable radius range of 213 
- 3O9i?0, under the assumption that the light variation is 
caused, now, only by the radius variation. In the studies of 
Haniff et al. (1995) and van Leeuwen ct al. (1997), i?*s of 
dozen of Mira variables arc presented, and most R*/Rq 
is larger than 300. The photosphcric radius of S Crt is 
close to the lower limits to those of Mira variables in the 
literature. Since the estimates of i?* and T BB are now 
available, the luminosity (L*) of S Crt is estimated to be 
2.29xl0 30 [W], resulting in the ratio L*/L of 5970. 

A similar analysis of data for a larger number of Mira 
and SRVs would be very useful for a better understanding 
of the difference and similarities between Mira variables 
and SRVs. 

5. Summary 

We have used VLBI monitoring observations with 
VERA to measure the absolute proper motions of a maser 
spot in S Crt with respect to the quasar J1147— 0724. We 
obtained an annual parallax of 2.33±0.13mas correspond- 
ing to a distance of 430^23 P c - 

Nearly two years astrometric observations revealed the 
internal motions of 26 maser spots for the first time. 
Taking into account the internal motions of the refer- 
ence maser spot, the absolute proper motion of S Crt sys- 
tem was obtained to be (—3.17 ±0.22 mas yr" 1 , — 5.41± 
0.22 mas yr _1 ), which is consistent with the result in the 
new HIPPARCOS catalog. 

The three-dimensional velocity field of the maser spots 
in S Crt was also detected, clearly outlining a bipolar out- 
flow with the major axis along northeast-southwest. The 
inclination of the axis with respect to the sky plane was 
~ 43°. The angular size of the maser distribution can 
be converted to the linear maser distribution radius of 
8.96 AU. 

Using the photosphcric temperature of 3097 K (Aringcr 
et al. 1999), the K band apparent magnitude mx =0.79, 
and the distance of 430 pc, the linear radius of S Crt 
photosphere (R*) was estimated to be 1.81 ± 0.14 x 
10 13 cm (260±20 Rq). This radius is comparable with the 
lower limits to those of Mira variables in literatures. 
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